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Abstract 

We present an investigation of the phase diagram of cholesteric liquid crystals within 
the framework of Landau - de Gennes theory. The free energy is modified to incorporate 
all three Frank elastic constants and to allow for a temperature dependent pitch in the 
cholesteric phase. It is found that the region of stability of the cubic blue phases depends 
significantly on the value of the clastic constants, being reduced when the bend elastic con- 
stant is larger than splay and when twist is smaller than the other two. Most dramatically 
we find a large increase in the region of stability of blue phase I, and a qualitative change 
in the phase diagram, in a system where the cholesteric phase displays helix inversion. 

PACS numbers: 61.30.Mp, 64.70.Md, 61.30.Dk 

1 Introduction 

Liquid crystals are anisotropic fluids typically composed of long, thin, rod-like molecules. They 
display long-range correlations in molecular orientation and show large length scale deforma- 
tions to even relatively weak external perturbations. As a result they may be well described 
at the continuum level by a vector field called the director ^ [2j , which represents the average 
local molecular alignment. 

A topic of long standing interest is the role of chirality in liquid crystals. The addition 
of small quantities of chiral dopant to a liquid crystal results in the appearance of a periodic 
structure in the molecular orientation with length scales typically in the optical range. The 
most common form for this periodic structure is a helical arrangement, known as the cholesteric 
phase, where the molecules display a natural twist along one direction. However, one of the 
most interesting features of chirality is that it also allows for more complicated structures. 
Although the cholesteric is always the thermodynamically stable phase at sufficiently low tem- 
peratures, it is found experimentally that upon cooling from the isotropic fluid, systems of high 
chirality display a series of first order phase transitions to brightly coloured 'blue phases' (BPs) 
before the cholesteric is reached. There are as many as three thermodynamically distinct blue 
phases in the absence of external fields, all of which appear only in a narrow temperature range 
(typically ~ IK) just below the clearing point. Two of them, BPI and BPII, display selective 
Bragg reflections in the visible range which can be indexed by cubic space groups, 0^~(/4i32) 
and 0^(^4232) respectively, while the third, BPIII, is amorphous. 

In the early and mid 1980's a series of important theoretical works showed how the cubic 
blue phases could be understood within both of the two principal theoretical continuum models 
of liquid crystals; the Oseen- Frank theory for the director field [31, and the Landau - de Gennes 
theory based on the Q-tensor order parameter |^ [3] (a review of both approaches is given in 
[S]). In both cases the picture which emerged was one of frustration between competing effects. 
From the director field point of view the locally prefered structure is one of ^double twisf, but 
this is found to be incompatible with global topological requirements. Therefore the cubic 
blue phases emerge as a regular array of double twist cylinders separated by a network of 
disclination lines. Complementing this the Landau - de Gennes theory picture of blue phases 
is that of a linear combination of biaxial helices chosen to optimise the competing bulk and 
gradient free energies. 
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The theoretical models were successful in accounting for the occurance, symmetry and 
general properties of the cubic blue phases, even assisting in the determination of specific 
space groups. However, the approximations necessary to make the analytic calculations feasible 
meant that quantitative comparison with experiment was largely not possible. Furthermore, 
the theories have a certain inflexibility since they retain a bare minimum of parameters and 
in this sense provide a 'one size fits all' description of the blue phases. This may have been 
adequate twenty years ago, but as the field has developed its shortcomings have become more 
apparent. One example of this is the insensitivity of the Landau - de Gennes theory of the 
blue phases to the value of the twist elastic constant. In contrast, the Frank director theory 
predicts that the twist elastic constant plays a central role [H] and furthermore, it has been 
observed experimentally that the stability of the cubic blue phases depends as much on the 
value of the twist elastic constant as it does on the chirality [Jj . 

More recently there has been considerable interest in trying to manipulate the properties of 
the cubic blue phases by adding a variety of chemical dopants [HI El 1^ 1^ > with the long term 
goal of utilising blue phases for device applications |121 1131 ITU I15j . One of the major obstacles 
to the use of blue phases in devices is their limited range of thermodynamic stability. However, 
Coles and Pivnenko have recently reported achieving a AOK range of stability for BPI in a bi- 
mesogenic compound doped with a chiral additive JHI- Although some questions remain as to 
whether their blue phase is thermodynamically stable or just metastable, it seems clear that 
the current theoretical picture of the blue phases is insufficient to account for all of the reported 
observations. Therefore, in this work, we ask whether extending the traditional Landau - de 
Gennes theory beyond a one elastic constant approximation will allow for a better comparison 
with experiments and perhaps aid the search for more practically useful blue phases. 

We first present a comparison of the Frank director and Landau - de Gennes theories, 
showing the relationships between them and discussing the limitations of both. Then we 
show how the Landau - de Gennes theory can be modified to account for three independent 
Frank elastic constants and a temperature dependent helical pitch in the cholesteric phase. 
The modified theory is treated in detail analytically for the cholesteric phase and it is shown 
that the theory allows for a description of the change in sense of the cholesteric helix upon 
decreasing temperature, which is observed in some systems P| I17L llSj . Next we present a 
numerical determination of the thermodynamic phase diagram of chiral liquid crystals for a 
range of physical parameters. To do this we introduce a new technique for determining the unit 
cell size of the blue phases, allowing for the first time, for a full minimisation of the free energy. 
Although the qualitative features remain unchanged when the elastic constants are varied, 
there is significant quantitative movement of the phase boundaries and the range of stability 
of the blue phases is brought to lower values of the chirality than has been reported previously. 
Our most striking result is that the range of stability of BPI is increased dramatically in 
systems where the cholesteric undergoes helical sense inversion. 

2 The Frank free energy and Landau - de Gennes theory 

The most widely used description of liquid crystals is the Frank director field theory ^ [^J. 
The local orientation of the liquid crystal is described by a directionless unit vector field n(r), 
called the director. It is found experimentally that the director satisfies an equivalence relation 
n ~ — n, so that its configuration space is the real projective plane, n S RP^. The static 
properties of liquid crystals are well described by director field configurations which minimise 
the Frank free energy subject to appropriate boundary conditions 
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where i^^i,i^22> and i^-ga 

are the Frank elastic constants, known as splay, twist and bend repec- 
tively. The parameter Qq is nonzero only in chiral liquid crystals where it determines the pitch 
of the cholesteric helix. 

In discussions of the cholesteric blue phases the Frank free energy is supplemented by a 
saddle-splay term ^ |2j 

(2.2) 

In ordinary liquid crystals this term can be safely treated as a total divergence and integrated 
to a surface term, which is then discarded. However, in systems possessing defects, such as 
the blue phases, variations in the magnitude of the order become as relevant as variations in 
the direction of the order. The value of the Frank elastic constants depends on the magnitude 
of the order (in a manner which can be determined using the Landau - de Gennes theory) so 
that in such cases the saddle-splay invariant may not be treated as a total divergence and can 
play an important role in the energetics of systems with defects. Indeed it has been shown 
that, within the director field framework, the saddle-splay term is entirely responsible for the 
thermodynamic stability of the blue phases El • 

There are two principal drawbacks of the director field theory; first, the magnitude of the 
order has to be put in by hand, and second, the effect of biaxiality is not included. Both of 
these shortcomings are most severe in the vicinity of defects where it is necessary to adopt 
a more sophisticated framework in order to give a complete description. This is a particular 
problem for theories of the blue phases since they contain a regular lattice of defects. An 
alternative option is to use a traceless symmetric second-rank tensor, Q, as a more general 
order parameter for liquid crystals. From a physical point of view it may be thought of as 
proportional to the anisotropic part of the magnetic susceptibility, or dielectric tensor. A 
phenomenological description of the energetics and phase transitions of liquid crystals can 
be provided by constructing a Landau theory using the Q-tensor 1 . This comprises a bulk 
contribution 

^buik Jd^r (atr(Q2) - 6tr(Q3) + c(tr(Q2))^) , (2.3) 

and a gradient contribution, accounting for the energy cost associated with distortions of the 
order, 

Fgrad^^ I ^3^ (^i^^(vxQ + 2goQ)' + ii2(V-Q)') . (2.4) 

The bulk free energy describes a first order phase transition between the ordered and disordered 
fluids. The parameters b and c are positive constants and a is a thermal scaling variable, which 
changes sign with decreasing temperature and has a linear dependence sufficiently close to the 
critical surface. The parameter qq defines the helical pitch. Since it is assumed that go is 
independent of temperature it follows that so too is the cholesteric pitch. It is also clear that 
this theory, taken to second order in the derivatives of Q, can only account for two independent 
Frank elastic constants; it has long been known that the constraint is that splay equals bend, 
Kf-^ = K^^. In order to remove this degeneracy, and allow for a temperature dependent helical 
pitch, it is necessary to retain higher order terms in the gradient free energy. 

The construction of invariants contributing to the gradient free energy at higher than 
quadratic order has been considered previously by a number of researchers |19( I2flj . Here, 
for completeness, we list the possible invariants which can be formed at cubic order in the 
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Q-tensor, and at most quadratic order in gradients. 



^abcQ adQ de^ bQ ce j 


(2.5a) 


(^abcQadQhe^ dQce j 


(2.5b) 


Qab^ aQbc^ dQcd , 


(2.5c) 


Qab^aQcd^bQcd , 


(2.5d) 


Qab^aQcd^cQbd , 


(2.5e) 


Qab^cQab^dQcd > 


(2.5f) 


Qab^cQoxNdQbd , 


(2.5g) 


Qab^cQad^cQbd , 


(2.5h) 


Qab^cQad^dQbc ■ 


(2.5i) 



We note that only six of the achiral invariants are linearly independent if, as is usual, we 
neglect total divergences, since we can construct the identity 

QcjN cQad^ dQbc + QcjN aQcd"^ cQbd - Q ah^ cQ ac^ dQ hd 

— Qab^aQbc^dQcd = ^ c(QabQad'^ dQbc — QabQ ac ^^dQbd) . (2.6) 

Since there are now many more Landau - de Gennes invariants than there are Prank elastic 
constants, it is instructive to see how the two descriptions are related, which may be done 
by letting the Q-tensor assume a uniaxial form, Qab = S{3nanh — Sab), and identifying the 
director field with the eigenvector of the Q-tensor corresponding to its maximal eigenvalue. 
A straightforward, but lengthy, calculation gives the following relations between the Q-tensor 
invariants and the director field: 



(VxQ)' 

(V-Q)^ 
Q V X Q 

^obcQadQde^bQ ce 
^obcQadQbe^dQ ce 
Qab^aQbc'^dQcd 

Qab^aQcd^bQcd 
Qab^aQcd^cQbd 

Qab^cQab^dQcd 
Qab^cQac^dQbd 

Qab^ cQad^ cQbd 
Qab^cQad^dQbc 



■ 18S2(Vanb)' + 5(Va5)' - QS\{V ani,){Vma) + (n X V X n)^) 

- 3(n -VSf - QS{VaS) (2n • Vua - riaV ■ n) , 

: 952((V • n)2 + (n X V X n)^) + 3(n -VSf + (V„5)^ 
+ 6S{VaS) (2n„V • n - n • Vn^) , 
: 95^ n • V X n , 
: 95^ n • V X n , 
: 95^ n • V X n , 

: 9S^(2(n X V X n)2 - (V • n)^) + S^Va(n ■ Vua - 2naV ■ n) 
+ 9S{n ■ VSf - S{VSf , 

: 185^ (3(n X V X n)2 - {Vamf) + 185(n • VSf - 6S(ySf , 
: 9^3 (2(n X V X n)2 - (V„n6)(Vbn„)) 

- S^Vai2n ■ Vua - UaV ■ n) + 95(n • VSf - S{VSf , 

: 185(n • VSf - QS{VSf - 6S^Va(n • Vn^ + n„V • n) , 
:9S3(2(V-n)2-(nx Vxn)2) 

- 2S^Va{^ ■ Vna + 4naV • n) + 9S(n • V^)^ - ^(V^)^ , 
QS{VSf + <9S''{Vanbf , 

: 9S3(2(Van6)(V6n„) - (n x V x n)^) 

- 253V„(4n • Vn„ + n„V • n) + 9S(n • V^)^ - S{VSf . 



(2.7) 

(2.8) 

(2.9) 
2.10) 
2.11) 

2.12) 

2.13) 

2.14) 

2.15) 

2.16) 

2.17) 

2.18) 
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Using the following identities for the director field 

i^amf = (V • n)2 + (n • V X n)2 + (n X V X n)^ + Va(n • Vn„ - n^V • n) , (2.19) 
(Vanb)(Vbn^) = (V • nf + Va(n • Vn„ - n„V • n) , (2.20) 

in Equations 1)2. 7() to ()2.18|) we can identify the Frank elastic constants in terms of the coeffi- 
cients of the Landau - de Gennes theory as 

2Kf^ = 9S\L2i + L22) + 953(-L33 - 2L34 - ^35 + 2L37 + L38 + 2L39) , (2.21) 

2Kl^ = 1852^21 + 953(-2L34 + L38) , (2.22) 

2Kl^ = 9S\L2i + L22) + 953(2L33 + 4L34 + 2L35 - L37 + - L39) , (2.23) 

4^12^0 = ^W^L2i + 953 (L31 + L32) . (2.24) 

It is not realistic to study the effect of all possible cubic terms in detail and consequently we 
shall restrict our attention to a subset. We choose to focus on three invariants, (|2.5aj) . ()2.5d|) 
and (l2:5hl) . These are chosen for a number of reasons; first, it is clear that one should retain 
at least one chiral and one achiral invariant. Since there does not appear to be any substantial 
difference between the two chiral invariants we feel confident that no qualitative changes would 
result from retaining both terms. It is more difficult to choose representatives of the achiral 
invariants. The invariant ()2.5h|) is retained because uniquely amongst the achiral invariants it 
contributes equally to all three Frank elastic constants. This leaves only a choice of invariant 
that will distinguish between the Frank splay and bend elastic constants. We choose to use 
(IT5d|) partly b ecause it has been used in previous work [211 122j . partly because it gives the 
largest distinction between Kf-^ and KJ^^ and partly because we found that it has the largest 
contribution to the energetics of a single double twist cylinder. 

To summarise, in the remainder of this paper we shall investigate the properties of the free 
energy 

F = V [d'rU tr(Q2) - ^/6tr(Q3) + (tr(Q2))' + ^ [(v x Q + Q)' + (V • Q)' 

^ ^ (2.25) 

+ g^Q2 . V X Q + ^QabVaQcdVbQcd + J^QabVcQadVcQbd] ) • 

We have made use of the well-known rescaling, Q (6/\/6c)Q, to reveal the redundancy of 
the parameters b and c in the bulk free energy and shifted to dimensionless variables. The 
change of variables is r ^ r/2qQ, F — > 6^/(288gQC^)F and we define 

L21 := 6cLi/b^ , L22 := 60^2/6^ , (2.26) 

K := 2qo^/L^ , r := 6ac/6^ . (2.27) 

We note that the parameter k, which is called the chirality, is the same as that defined by 
Grebel et. al. j^, but that our reduced temperature r differs from their definition; the two 
being related according to ^ghs = t + k^, where tons is the reduced temperature used in |^. 



3 Cholesteric Phase 

We first investigate how including cubic order invariants in the gradient free energy. Equation 
(|2.25|) . modifies the usual theory of the cholesteric helix. It will be assumed that the helical 
order parameter still takes the form dictated by the quadratic theory, namely 

Q = ^^(0 -1 |-h%fo cos(A;a;) sin(A:x) I . (3.1) 




V2 



i^O sin(A:x) 
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With this form of Q-tensor the expression for the free energy becomes [SHj 



+ xQ2[(l- Wfc'-2(l-aQ,.)fc] , (3.2) 

where, a := — -L3i/(-v/6-^^2i) and /3 := (2L34 — L38)/(\/6-^^2i)- Minimising the free energy leads 
to the set of Euler-Lagrange equations 

^ = i-^Qft ' (3.3a) 

\(t + /^2)g^ + 3(Q2 _ g2) ^ 4g^(g2 ^ g2) ^ ^g2 [_^^2 ^ 2aA;] = , (3.3b) 

\(t + Ac2)Q2 - eg^^Qs + 4Q2(Q' + Q2) + ^ Q2 [(1 - mh)k^ - 2(1 - aQ/,)A;] = . (3.3c) 
Using these, and assuming 7^ as expected for a first order transition, we readily obtain 



Ql = \ (l2Q. + ^^^^^^ - (r + K^)) - Q 



2 



(3.4) 



+ -^Q.^^^^^ = (e - * [2« - /3lEfg^] ^EfS^) . (3.5) 

In general these equations must be solved numerically; however, they can be easily solved when 
a = /3 = 0, and we recover the usual results ^ 101 EH 



Ql = '-^fQh - i . (3.6b) 

The main feature which arises from the inclusion of higher order terms is that the heli- 
cal wavevector now depends on the temperature. As Equation 1)3. 3a() shows, the wavevector 
depends on the amplitude of the homogeneous part of the order parameter, and since the 
order parameter depends on the temperature, it follows that so does the cholesteric pitch. It 
should be noted that, even with such a simple modification of the free energy, we can already 
account for the experimentally observed inversion of the helical sense with decreasing temper- 
ature found in some cholesterics ^ I17| I18j . According to the present theory the cholesteric 
wavevector changes sign when = 1/a. When this value is substituted into Equations 1)3.4(1 
and ()3.5|1 we obtain the following expression for the inversion temperature 

THi = ^(3a-8) . (3.7) 
At the inversion itself we find Q2 = V^Qh and 



Q = — 7= 2 , (3. 




which corresponds as expected to a uniaxial nematic with director n = (0,1,0). Other ori- 
entations of the nematic within the yz-plane may be obtained by adding a constant phase 
shift to the sin and cos terms in Q in Equation ()3.H) . It is clear that within this theory, as 
observed experimentally, the helical inversion is a smooth transition which does not involve 
any discontinuities in physical quantities, i.e., it is not a phase transition. 
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(a) 



(b) 



Figure 1: (a) Isotropic cholesteric transition temperature as a function of chirality for a range 
of parameter values: (i) a — (3 ^ 0, (ii) a = 2 , /? = 0, (iii) a = 2 , /3 = 1, (iv) a ~ i , P ~ 0. 
(b) Temperature dependence of the cholesteric wavevector for the same set of parameters as 
in (a) (all cases are for k = 0.05). 

The isotropic-cholesteric transition temperature is obtained by supplementing Equations 
H3.4|) and (|3.5|) with the condition that the free energy (Equation ()3.2|) ') be zero. A short 
calculation leads to the following equation for Qh at the transition temperature: 



1-aC 



2 1-fSQh 



24Qh 



2a -pj 



l-aC 



-PQh 



1-aC 



2a -Pj 



l-aC 



Again when a = /3 = we recover the usual results |1J IHl 123! 



k2n3/2 



. (3.9) 

(3.10a) 
(3.10b) 



where tic is the temperature at which the isotropic-cholesteric transition occurs. A situation 
of some experimental interest occurs when the cholesteric undergoes helix inversion at the 
isotropic-cholesteric transition temperature ^7j- This case may be solved exactly by setting 
Qh = 1/a in Equation (|3.9j) . We find a = 4 and 



''"HI = nc 



1 



(3.11) 



which is precisely the transition temperature for a nematic (with non-zero k). This particular 
case provides a guide to the range of values that a can be expected to take in experimental 
systems. 



In the general case the equation can be solved numerically. Figure 1(a) shows the isotropic- 
cholesteric transition temperature as a function of chirality for a selection of values of a and 
/3. The general trend can be easily understood on the basis of the effect that the parameters 
have on the helical pitch. Increasing a increases the pitch and shifts the liquid crystal towards 
nematic behaviour, while increasing /3 decreases the pitch and hence shifts away from the 
nematic state. The temperature dependence of the helical wavevector is shown for the same 



set of values of a and (3 in Figure 1(b) It can be seen that the wavevector varies approximately 
linearly with temperature except in the immediate vicinity of the transition temperature. It 
may at first seem unusual that increasing [5 should lead to a larger temperature dependence, 
however this can be understood by noting that the inversion temperature depends only on a 
and that increasing (3 leads to a larger wavevector at the transition temperature. 

Finally, it is known that if the chirality is sufficiently large then the isotropic-cholesteric 
transition is second order instead of first. We can calculate the maximum value of k compatible 
with a first order transition by setting Qh = in Equation (jSIHI)- The resulting equation is 
solved to obtain 

_ 2 



-1 + + 3(2a - /?) 



(3.12) 



"'max — 2a-l3 

In what follows it will be seen that this limit is never reached due to the intervention of blue 
phases. 



4 Numerical Method 

Our next aim in this paper is to construct the phase diagram resulting from a minimisation 
of the modified Landau - de Gennes free energy. Equation (|2.25|) . taking into account the two 
cubic blue phases. Although a certain amount on the blue phases can be done analytically, 
this approach involves adopting an approximate form for the order parameter and only gives a 
constrained minimisation. Furthermore, althougth the analytic theory correctly identifies the 
two structures observed in experiments, it does not reproduce the correct order of appearance 
of the two phases at low chiralities. A full minimisation can be achieved numerically and 
has been described recently in Reference ,26^ for the one elastic constant approximation. The 
Q-tensor is relaxed towards the minimum of the free energy according to the equation |27j 

^ = THa, , (4.1) 
where F is a collective rotational diffusion constant and the molecular field is given by 

Since we are entirely concerned with static equilibrium configurations we have neglected the 
coupling to fluid flow, both in the above equations and in our numerical simulations. The 
equations are solved using a three dimensional lattice Boltzmann algorithm, the details of 
which have been given in Reference j28j . 

To study the different phases, cholesteric, BPI or BPII, it is necessary to implement appro- 
priate initial conditions for the simulation. The Q-tensor is initialised using analytic expres- 
sions appropriate to the high chirality limit which act to define the symmetry of the chosen 
phase. Under subsequent numerical evolution according to Equation (|4.1|) the system relaxes 
to that structure of the same symmetry which minimises the free energy. We are therefore 
able to obtain, for any value of the parameters, local minima of the free energy corresponding 
to each of the cholesteric and blue phases. The global free energy minimum was taken to be 
the smallest of these calculated local minima. 

We have seen in Section |21 that the inclusion of cubic invariants in the gradient free energy 
leads to a temperature dependent helical pitch in the cholesteric phase. Equation 1)3. 3a() . For 
the blue phases as well, the unit cell size is temperature dependent, so that to achieve a full 
minimisation of the free energy it is necessary to set the correct unit cell size in the simulation. 
This unit cell size is not known a priori, but rather depends on the magnitude of the order 
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-4 -3 -2-1 1 -4 -3 -2-1 1 



(a) (b) 

Figure 2: Comparison of lattice Boltzmann (o) and analytic (solid line) solutions for the 
helical pitch in a cholesteric displaying helical sense inversion, (a) Cholesteric wavevector 
and (b) pitch against reduced temperature, r. 

parameter, a quantity which is only determined by the numerical minimisation. Therefore 
we must introduce a means of determining, and setting, the unit cell size as the Q-tensor 
evolves during the simulation. We can account for a change in unit cell size by rescaling the 
gradient contributions to the free energy and molecular field. This is accomplished in practice 
by changing the elastic constants as follows 

qo = <itlr , (4.3a) 
L2a = X , (4.3b) 
= 4f X r2 , (4.3c) 

where a = 1,2, b = 1,...,9, a superscript 'init' denotes the initial value of a simulation 
parameter and r is the appropriate rescaling factor, which in previous analytic and numerical 
j26j work was referred to as the 'redshift'. One problem with the analytic theories is that the 
value of the redshift is not determined exactly, but only for the approximate form of the Q- 
tensor that is assumed. Similarly, in previous analytic work, the redshift was assumed to 
take the value suggested by the approximate analytic calculations. The exact redshift for the 
cholesteric phase could be calculated by obtaining the numerical value of Qh and using Equation 
(|3.3a|) . However, a similar approach is not available for the blue phases and consequently it is 
more useful (and easier) to calculate it using the free energy as follows: since the free energy 
is quadratic in gradients, it may be written formally in k-space as 

F = ak'^ + bk + c, (4.4) 

where the coefficients a, h and c depend on the Q-tensor, but not on k. The optimum wavevector 
is given hy k = — 6/2a, and since the coefficients a and h are determined by the simulation 
every timestep it is straightforward to use these values to determine the exact value for the 
redshift. 

In order to verify that the procedure was working successfully, and to check the level of 
accuracy that could be obtained, we used the lattice Boltzmann algorithm to calculate the 
wavevector of a cholesteric undergoing helix inversion and compared it to the theory described 
in SectionlSI Simulation parameters were chosen to set a = 2.0 and k ~ 0.096. The results are 
shown in Figure |21 We have plotted both the helical wavevector, which is the relevant quantity 
theoretically, and the pitch, since this is more frequently given in experimental work. As can 
be seen the agreement is excellent even very close to the inversion point. 
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1.2 1.4 1.6 



(a) 



(b) 



Figure 3: Numerical phase diagrams in the one elastic constant approximation, r, the reduced 
temperature, and k, the chirality, are defined by Equations (|2.27|l . (a) Optimisation of the 
unit cell size, (b) Phase diagram with fixed unit cell size, after |26|. 

4.1 Results 

In this section we describe the numerical phase diagram obtained using the lattice Boltzmann 



algorithm outlined above. Figure 3(a) shows the phase diagram obtained for the one elastic 
constant approximation with numerical optimisation of the unit cell size. For comparison, the 
phase diagram for a fixed unit cell size determined in Reference [SHj is reproduced in Figure 
|3(b)| We note that optimisation of the unit cell size has extended the range of stability of BPI, 
both relative to the cholesteric phase and relative to BPII. The movement of the cholesteric 
phase boundary is quite significant, with the triple point moving to lower chirality by about 
20%. This is due to the optimum redshift taking a lower value at these chiralities than was 
assumed previously. Based on analytic calculations of Grebel et. al. [S] the redshift was 
assigned the value 0.79 in j26j . However, we find a much lower value, with average 0.68 at the 
cholesteric-BPI phase boundary (the value is roughly independent of temperature, except very 
close to the isotropic transition temperature). In contrast, at the BPI-BPII phase boundary 
we obtain an average redshift of 0.77 for BPI and 0.86 for BPII. These values are in better 
agreement with the analytic results of Grebel et. al, primarily because the chirality at this 
phase boundary is much closer to the values of chirality where the analytic calculations predict 



the blue phases are stable. In what follows we will use Figure 3(a) as a reference phase diagram 
relative to which the effect of varying the Landau - de Gennes parameters can by measured. 

We next constructed phase diagrams for different ratios of the elastic constants. Here 
one expects that the qualitative features of the phase diagram will be retained, but it is 
nonetheless of interest to determine how large a quantitative shift can be obtained. Figure 
|1] shows the phase diagrams obtained upon separately varying the twist and bend elastic 
constants. To investigate the effect of the bend elastic constant we chose parameter values 
L21 = L22 = Lm = 0.02, L38 = which corresponds to a ratio of bend to splay of about 
1.75, while splay and twist remain degenerate. The stability of BPI is seen to decrease quite 
significantly relative to the cholesteric phase while at the same time there is a small increase 
in stability over BPII. There is only a minor shift in the cholesteric-BPI phase boundary at 
the transition temperature, however, as the temperature decreases the shift becomes larger. 
For example at a reduced temperature of r = —2 the phase boundary occurs at a chirality of 
K ~ 1, representing a shift to higher chiralities of almost 70% as compared to the one elastic 
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Figure 4: Phase diagrams obtained for different values of the elastic constants, r, the reduced 
temperature, and k, the chirality, are defined by Equations (|2.27() . (a) Kf — ifj — 0.5K. 



3 ' 



(b) if^ =X| =1.5if|. 



constant case. The value of the redshift is markedly increased, with the average value for 
BPI at the cholesteric boundary being 0.81. However, it was shown in Section |21 that the 
cholesteric wavevector depends on the value of L34 and hence has also increased. This raw 
value of the redshift no longer represents the most directly accessible quantity and a more 
relevant figure is the ratio of the BPI redshift to that of the cholesteric, since this can be 
measured experimentally by means of the discontinuity of the back-scattered Bragg peak. For 
this quantity we obtain an average of 0.67 and note the strong similarity of this value with 
that obtained in the one elastic constant approximation. 

At the BPI-BPII phase boundary the redshift is 0.87 and 0.97 for BPI and BPII respectively. 
Again, it is the ratio of these values which is more directly relevant to experiment. The ratio 
of the unit cell size to that of is 0.90, and again we note a strong similarity with the 
ratio obtained from the one elastic constant approximation, 0.89. 

The value of the twist elastic constant is controlled by the Landau - de Gennes parameter 
L22- In most liquid crystals the twist elastic constant is smaller than either splay or bend. In 
order to match this, we constructed the phase diagram for parameter values L21 = 0.02, L22 = 
0.04, L34 = L38 = 0, which is shown in Figure p(b)[ This choice of parameters resulted in a ratio 
of splay to twist of about 1.5, while splay and bend remained degenerate. Again we observe that 
the stability of BPI is reduced relative to the cholesteric phase by an amount similar to that 
seen by varying the bend elastic constant. The values of the redshift for both BPI and BPII 
are only very slightly increased relative to their values in the one elastic constant limit, while 
the cholesteric wavevector is insensitive to the value of L22- This reveals an intriguing feature, 
that while the phase boundaries and absolute values of the redshift can vary appreciably, the 
ratios of the redshift for the different phases, and hence the discontinuities in back-scattered 
Bragg peaks, are essentially independent of the values of the elastic constants. 

Finally, we investigated the effect of the chiral cubic invariant on the blue phases. We 
chose to use L21 = L22 = = 0.02, L38 = 0.04, a = 2.0. (For the cholesteric phase this sets 
P = and gives a ratio of bend to splay of about 1.6.) The phase diagram is shown in Figure 
El The value of a is such that the cholesteric undergoes helical sense inversion at a reduced 
temperature of about r = —2. What is remarkable is the dramatic increase in stability of BPI 
relative to the cholesteric phase. The region of stability has been increased down to chiralities 
as low as K = 0.07 and at such low chiralities the phase boundary is essentially independent 
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Figure 5: (a) The numerical phase diagram obtained with the chiral invariant, H2.5a|l . added 
to the free energy. The magnitude of this term was chosen so as to produce heHcal sense 
inversion in the cholesteric phase at a temperature not far below the isotropic transition 
temperature, r, the reduced temperature, and k, the chirality, are defined by Equations 
(|2.27ll . (b) An enlargement of the region near the isotropic transition temperature. Note the 
reversal in the order of appearance of BPI and BPII as a function of chirality. 

of K for all r. In addition, we find a very small region of stability for BPII located close to 
the isotropic transition. As shown in greater detail in Figure 5(b) this occurs in a narrow 
temperature interval at chiralities lower than those for which BPI is stable, representing a 
reversal of the order of appearance of the two blue phases. In the region where BPII is stable 
we find a redshift of 0.36. In contrast to the situation without the chiral invariant, this is 
smaller than the BPI redshift, which takes the value 0.43. Again, the absolute values of the 
blue phase redshifts are not as relevant as their ratios to that of the cholesteric phase. In this 
case we find the cholesteric redshift is 0.49 at the isotropic transition temperature, giving ratios 
of 0.88 for BPI and 0.74 for BPII, both of which are significantly different to those obtained 
in the one elastic constant limit. 

Since BPI is now stable over a much larger temperature range it displays a significant 
variation in unit cell size as the temperature is lowered. As an illustration of this, the BPI 
redshift is 0.19 at a reduced temperature of r ~ —5, corresponding to more than a two-fold 
increase in the lattice parameter. A plot of the temperature dependence of the BPI redshift 
is shown in Figure El For comparison the helical wavevector of the cholesteric phase is also 
plotted on the same graph. 



5 Conclusion 

We have investigated numerically the phase diagram of the cholesteric blue phases for a range 
of parameter values within the framework of a modified Landau - de Gennes theory. The 
traditional Landau - de Gennes theory has long been known to only accommodate two inde- 
pendent Frank elastic constants and to have a temperature independent cholesteric pitch. Both 
of these shortcomings were overcome by retaining terms of cubic order in the Q-tensor in the 
expansion for the gradient free energy. Since the new terms were added specifically to remove 
the degeneracy between splay and bend, and to give a temperature dependence to the helical 
pitch they possess a clear and simple physical interpretation. In particular the value of the 
parameter a which controls the strength of the chiral cubic invariant should be relatively easy 



12 





(a) (b) 

Figure 6: Lattice Boltzmann results for the temperature dependence of the (a) wavevector 
and (b) pitch of the unit cell of BPI along the cholesteric-BPI phase boundary at k w 0.07 
(o) . For comparison the cholesteric wavevector and pitch are also given: lattice Boltzmann 
results (o) , results from the calculation presented in Section |31 (solid line). 

to estimate on the basis of Equation ()3.7() for the helix inversion temperature. The magnitude 
of the elastic constants for the achiral cubic invariants is more difficult to ascertain. Ahhough 
this may be estimated from the ratios of the Frank elastic constants it is clear that because 
there are more Landau - de Gennes elastic constants than Frank, the latter are insufficient 
to uniquely determine the former (an estimation of the magnitude of the cubic Landau - de 
Gennes elastic constants was given in jl9jl. 

The modified Landau - de Gennes theory that we have investigated provides a phenomeno- 
logical description of helical sense inversion in the cholesteric phase. The inversion arises as a 
natural consequence of the presence of including higher order chiral invariants in the gradient 
free energy. Although such higher order contributions are usually neglected since they are 
deemed small compared to the terms already retained, the fact that helix inversion is observed 
experimentally demonstrates that these terms can play a significant role. Mathematically, we 
comment that for systems undergoing first order phase transitions (such as liquid crystals), 
although the order parameter is small it is not infinitesimal, and therefore the relative mag- 
nitude of a given term in the Landau expansion depends not only on the power of the order 
parameter but also on the size of any numerical coefficient premultiplying it. 

The primary aim of this paper was to investigate how much the properties of the blue 
phases could be changed within the framework of Landau - de Gennes theory. In this regard 
we have shown that the retention of cubic order terms in the gradient free energy can lead 
to considerable changes in the size of the blue phase unit cell and in their phase diagram. 
Most dramatic amongst the results is the increase in stability of blue phase I obtained in 
systems where the cholesteric undergoes helical sense inversion. Again we comment that this 
significant, qualitative change in the phase diagram arises from retaining cubic order terms and 
demonstrates that these give rise to more than just small changes in the physical properties. 

It is of interest to consider whether the mechanism considered here is a candidate to account 
for the increased range of stability in blue phase I recently reported by Coles and Pivnenko jl6j . 
It seems not, as apart from the large temperature range most features of their blue phase differ 
from those obtained for the choice of parameters we made here. For example, the numerics 
show an increase in the size of the BPI unit cell with decreasing temperature, while in the 
experimental system the unit cell size shows a small decrease. Also, numerically we find that 
BPII has a larger unit cell than BPI at the transition between the two, in contrast to what 
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is found experimentally. However, we remark that the parameter space in the Landau - de 
Gennes theory is large enough that these discrepancies could well be resolved by a different 
choice of parameters. 
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